REDUCTION OF THE N-COMPONENT SCALAR MODEL AT TWO 

LOOP LEVEL 



Antal Jakovac 
Department of Atomic Physics 
Eotvos University 
Budapest, Hungary 
(February 7, 2008) 

Dimensional reduction of high temperature field theories improves IR features of 
their perturbative treatment. A crucial question is, what three-dimensional theory 
, is representing the full system the most faithful way. Careful investigation of the 

ON ■ induced 3-dimensional counterterm structure of the finite temperature 4D O(N) 

symmetric scalar theory at 2-loop level leads to proposing the presence of non-local 
operators in the effective theory. On scales beyond 0(T~ l ), the scaling behavior 
("^ ■ of the couplings, consequently, deviates from the usual three-dimensional scaling 

\ characteristic for superrenormalizable theories. 
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1 I. INTRODUCTION 



Numerical investigation of the finite temperature electroweak phase transition is unavoidable 
because in the symmetric phase no improved perturbation method is known to work reliably. The 
most adequate approach is to concentrate on the degrees of freedom left "massless", e.g. the static 
I/"") . Matsubara modes of the magnetic gauge and Higgs fluctuations (T) . 

The effective reduced theory usually is arrived at after simple f-loop integration over the non- 
static and the screened static modes ^,|| , with the restriction that only those terms are retained in 
the effective action which are renormalizable in 4 dimensions. These theories being superrenormal- 
izable in 3 dimensions one can work out their exact divergence structure and relate the physical, 
temperature dependent mass to the bare parameters of the 3 dimensional theory . The general 
strategy behind this procedure is the matching of some important perturbatively computable quan- 
tities calculated both in the full finite temperature theory and in the effective superrenormalisable 
three-dimensional theory fi. The difficult and most interesting question is how accurate is this 
3 dimensional superrenormalizable representation of the theory. On one hand one should control 
higher dimensional operators with couplings inversely proportional to some power of T |}| . On the 
other hand on very general grounds intrinsic non-local behaviour of range T~ 1 is also expected to 
appear. 

The last point can be illustrated, for instance, on the example of the so-called "sunset" diagram 
contributing to the propagation of a static mode. The contributions are naturally divided into 
three groups: 




--0 



(1) 



(the solid lines represent non-static, the dashed ones the static propagators). The first class is 
accounted for when one integrates over the non-static modes with two-loop accuracy. The second 
class is considered in the course of the solution of the static effective model. The third ("mixed") 
class, however would be left out from the two-step solution. 

Our proposition is to incorporate the mixed part of the "sunset" contribution into the effective 
model by defining a new kind of vertex, namely 
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This vertex is clearly nonlocal, since the lines compressed into it have a typical range C(T _1 ). 
Similar conclusion has been emphasized recently also by [Q. 

Explicit estimate of the strength of the non-local corrections, to our knowledge will appear for the 
first time in the present paper. In this first attempt we limit our calculations to the N-componcnt 
scalar model (for N = 4 the Higgs-sector of the standard model) in the hope, that the main lesson 
remains true for the full electroweak theory || . 

We did our computations using momentum cut-off regularization, that is we regularized the prop- 
agator to be zero, when its momentum exceeds the cut-off. The divergent counterterm structure 
obtained for the effective model by the 2-loop reduction and that of the 3 dimensional superrenor- 
malisable theory computed with cut-off regularization both display extra power divergences which 
would be absent in dimensional regularization |^,| 10 1. The need for extra pieces in the effective 



action becomes obvious by the mismatch discovered just between those quantities. The correct 
choice of the new terms will be undisputable if with their contributions the balance between three 
different type of power- and logarithmic divergences of the effective theory and those induced by 
the reduction will be reestablished. 

Our investigation starts with a detailed presentation of the counterterm structure of the 3 dimen- 
sional O(N) symmetric model of the static modes induced upon 2-loop integration over nonstatic 
degrees of freedom (Section 2). This is then compared to the divergences of the effective potential 
calculated from the 3 dimensional local superrenormalizable O(N) model at the same 2-loop level 
(Section 3). A discrepancy between the results will be discovered and the origin of the mismatch 
will be located in Section 4. The extra divergencies needed for the consistency between the effective 
and the original theory come from intrinsically non-local operators in 3 dimensions with a charac- 
teristic non-locality range T _1 , as qualitatively explained above. In the Conclusion (Section 5) we 
shall outline the matching strategy for replacing the effective non-local theory by a local theory 
with nonremovable cut-off aT^ 1 (a — 1 — 5). The interpretation of the computer simulation of a 
corresponding lattice system will also be shortly discussed. Also we shall touch upon the extension 
of the present work to gauge theories. 



II. REDUCTION OF THE 4 DIMENSIONAL O(N) MODEL AT TWO LOOP LEVEL 

The model is described by the following Hamiltonian: 



N 1 \ ( N \ 2 
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(3) 



where H c t represents the countertems. First, we integrate out the nonstatic Matsubara modes, 
and find the potential energy density for the reduced 3 dimensional theory. For this we shift the 
fields ip by a constant <&o, what stands for the static part. Because of the O(N) symmetry there 
exists a coordinate system, where $o points along the N-th axis. The Hamiltonian breaks up into 
(^-independent, and quadratic pieces, and a higher power ^-dependent term: 



H = H a +H 2 +Hi, (4) 



where 
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We compute the potential energy density of the effective theory / in the following way: 

/y = J ff + ^TrIogA-(e-^-l> c , 



(6) 



where the index c refers to connected graphs, V is the volume of the system. 

For two loop results the interaction part of the exponential must be expanded to 0(X 2 ), and 
only <j> -dependent terms should be retained: 

/N-l \ 2 AT— 1 

( E^ 2 )c+ 2 (^E^ 2 )c+(^)c 



24 



1=1 

N 



(7) 



i=l i=l 

This result can be rewritten in terms of Feynman-graphs. We introduce the following notations 
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J(m) := — Tr log A, 
K(m 1 ,m 2 ,m 3 ) := 




(8) 



where mi, m2, m.3 are the masses of the propagators on the three internal lines of the diagram. It 
is clear, that 




dl(m) 
dm 2 



=:I'(m). 



(9) 



We introduce the notations 



m% := m 2 + ^$ 2 , 
m 2 H := to 2 + ^$ 2 , 



(10) 



and get the following result: 



*t + \(N l)I(m G ) + \l{m H ) + \{5m 2 + y )(iV - l)J'(m G ) 



+ I ((5m 2 + ^)j'( mjr ) + A [(TV 2 _ l)/'(m G ) 2 + 2(7V - l)/'(m G )/'(m H ) + 3/'(m ff ) 2 
A 2 

~ 36 $ ° [3^( m ff : rn H ,m H ) + (N - l)K(m H , m G , m G )] 



(11) 



The functions / and K can be expanded with respect to the masses, if no IR divergencies arise. 
This happens in our case, where the IR sensitive part is substracted because no static mode is 
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allowed to propagate on the internal lines. Therefore, one is allowed to expand I and K in powers 
of to 2 /T 2 . If we don't want to keep the operators suppressed by some inverse power of T, we can 
truncate the expansion after the first few terms (high temperature expansion): 

7(m) = I + hm 2 + I 2 m 4 + 
I'{m) =h+ 2I 2 m 2 + 3/3TO 4 + . . . , 

m? + to 2 . + to? 

K(m 1 ,m 2 ,m 3 ) = K + K 1 ^ ^ 3 + ... (12) 

(the function K is symmetric in the three m's). Relying on the divergence structure of the zero 
temperature model, and on dimensional analysis, the expected form of the coefficients in ( |l2] ) is: 

rl a 2 1 rl^A 1 T0rp"2 



h = If A 2 + J-fTA + T{T l 
T, 
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K Q = K 2 A 2 + KqTA + K l oa T 2 log f - ) + K°T 2 , 
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X: = iff °» (fog - j + log I —J + K». (13) 

For future convenience let us introduce N — (N + 2)/3. At one loop level, if we choose the 
counterterms (in general renormalization scheme) as 

Sm 2 = Am 2 ~ NX Q/ 2 A 2 + m 2 I l 2 ° 9 log (± 

5Ai = AAi - (N + 2)X 2 I l ° 9 log f - J , (14) 

where the terms A are finite (and specify the renormalisation scheme), we get for the effective 
potential: 



f = m + Am i + NX ( § 7 i T + % m ) + N \m ±2 log j + Y A^TA 



^ (A + AAi + (N + 2)X 2 I° 2 +(N + 2)X 2 I l ° 9 log (|) ) + . . . (15) 
At two loop level, with the same procedure we choose the counterterms as 
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6m 2 = Am 2 2 NA 2 l^AXJ 2 - \ K 2 X 2 - *±1]*%°\* log 

+N ((A> + 2)(4° 9 ) 2 + jUf A 2 m 2 (log^^ 
+N ( 2 (/^) 2 + ijjf*) A 2 m 2 log log (|) 
-A> (Am^A + AAi/^TO 2 - 2I l 2 ° 9 1«X 2 m 2 - U<[ 09 X 2 m 2 ^ log ^ 
SA 2 = AA 2 + 20 + 22iV + 3A> 2 (4og)2 + 5iV + 4^ os A ,3 A AX 2 
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Using also the expression of the one loop counterterms (|TJ) we obtain the 2-loop result valid in any 
general renormalization scheme. Some important systematics can be discovered in this expression 
(see the comments below eq.(p2|)), if one introduces the following shorthand notation: 



(17) 



A = A + AAi + 7V/ 2 °A 2 + (N + 2)I l 2 ° 9 X 2 log (£) 



Then the two-loop effective potential, renormalised from the point of view of four-dimensional 
ultraviolet behavior can be written with 0(A 3 ) accuracy as 



K + Am^)(l + 7VA(/ 2 +4 O9 log(^ 



2#*« + l -K l A log (£) - Ul^f + Ik*°>) (log £) 



+i VT 2 A ( -J° - -K° X - NX(I°? 2 ° a + -iC) log ( - 
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If we chppse = T, to avoid large logarithms, and a special scheme, eg. the following (Weinberg- 
type) one a 



d 2 f 
9$ 2 

d 4 f 



:= to 2 , 



A, 



(19) 



M=T=0,* =0 

we can obtain a considerably simpler result. At one loop level 



/ = ( «<" 



-NXI?T 2 + -NXllTA 
2 1 2 1 



24 



rX^n 



(20) 



which corresponds to the following choice of the finite parts: 



1 The correct way would be to choose fj, to be equal to the scale where the parameters A and m are fixed 
(eg. the mass of the Higgs boson in the present case). The large logarithms can be summed with the help 
of RG. This would give the above result with logarithmically T-dependent A and m. 
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Ami = -N\I$m 2 , 

AAi = -(7V + 2)A 2 / 2 °. (21) 

One might notice, that at 1-loop the coefficients of T 2 and of AT do not depend on the details, 
how the regularising cut-off is imposed. 

At two loop level (without explicitly giving the complicated expressions for Am\ and AA2) 

' "'- 4 w " 1 1 \ ■ 1 • • ■ (| --u v • , ■■• ( ...„ j ■ v T 



+NAT [\XI\ A 2 + l -Kl + log (± 



Y A X ^ ( 22 ) 



It can be seen from formula (18), that the coefficients of the terms 3> 2 T 2 and $ 2 AT do not depend 
on the renormalisation scheme and scale (/i). That means, that essentially this form would have 
arrived also at any other renormalisation scheme. On the other hand the coefficient of $ 2 T 2 and 
<& 2 AT at C(A 2 ) depend on the specific implementation of the cut-off (regularisation dependence), 
through the coefficients K\ and K®. Therefore for computing physical quantities it is necessary 
to apply a common unique regularization procedure, otherwise the 3D linear divergence would 
not be cancelled and also some C(A 2 ) finite contributions would be unreliable. The coefficients of 
other induced divergencies of the type T 2 log(A/T), ATlog(A/T) do not depend on the specific 
implementation of the cut-off in multiloop integrals. 

The only remaining task is the computation of the constants defined in ([l3"|). The Ii quantities 
come from the tadpole graph: 

h + 2I 2 M 2 + 3/3M 4 + ...=TJ2 [ ffv 2 ^ \—T72 - (23) 

$£qJ ( 2?r ) <+P +M 2 

The Ki constants are determined by the setting sun diagram: 

K + K\M + K 2 M + . . . —T J_^J (27r)9 { «l + El) {< + Et) { u> 2 + Eiy ™ 

where E 2 = p 2 + M 2 . The result of the integrations (see Appendix A and B) are: 
A 2 AT T 2 



h 



12 = —nr^ l °g T7Z2 (! + lo s( 27r ) " 7b), 



8tt 2 2tt 2 12 : 

i6^ los lr; Ki,T- 



h = 



C(3) 



192tt 4 T 2 ' 



K = 0.0001041333 A 2 - 0.0029850437 AT + -^T 2 log ( - ) - 0.0152887686 T 2 , 



32tt 2 "°\T 

Ki = -jjj^l ( lo S + °- 001087871 lo S (f) + const ( 25 ) 

The coefficients written in decimal form are results of numerical integrations, other coefficients were 
calculated analytically. The former depend on the special implementation of the cut-off procedure 
the latter do not. 

With these values we obtain the following result for the 3 dimensional effective potential in 
Weinberg-type renormalization scheme: 

A 2 A 2 /A' 

- -0.0013551443 A 2 log - 



kYA.T ( - A + 0.0012227544 A 2 - A_ iog (± 



±K (26) 
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The one loop result is well known [|2|j2|], the two loop result is new. From this expression one can 
read off the induced 3 dimensional mass counterterm: 

Sml d = -N _ 0.048277409 A) AT + ^ AT log (A) + ^T 2 log (j^) (27) 

No new divergence related to A appears, supporting the superrenormalizable nature of the 3 di- 
mensional theory. We have to compare this induced counterterm to the counterterm required for 
the 3 dimensional theory of the same form as (||) expected on the basis of superrenormalizability. 
The eventual cancellation is a precondition for the consistent representation of the 4D finite T 
theory by this simplest minded 3D form. But as we shall see this hope doesn't come true. 



III. 2-LOOP EFFECTIVE POTENTIAL OF THE 3D HIGGS-MODEL 



The result (26) of the previous section suggests a 3D effective representation of the original theory 
which could be of the form (|^). This is again a local scalar O(N) model, just now with temperature 
dependent mass parameter. The renormalization of the theory, however, on the 4D level has 
been accomplished already. The counterterms of the "embedded" 3D theory were generated by 
the reduction. A crucial check of consistency on the suggested model is the cancellation of its 
divergencies against the induced counterterms. 

In this section we investigate the effective potential of the three dimensional theory defined by 
the action (||) at two loop level, when the bare mass is temperature dependent, as described by 
(|26]). The results are known for even more complex theories so this section will be just a short 
summary of well-known facts. 

In three dimensions the theory is superrenormalizable, we expect divergent counterterm only to 
the mass parameter, and only up to two loop level. Otherwise the graphs are the same as in the 
4 dimensional case (see (pi)), but the integrals are 3 dimensional. The result of the integration is 
(cf. ©): 

/ = \{™ 2 + 5m 2 )<S>l + ±$t + \{N - l)I(m G ) + ±I(m H ) + l -5m 2 {N - l)I'(m G ) 

+ ^Sm 2 I\m H ) + A [(iV 2 - \)I'{m G f + 2(A^ - l)/'(m G )/'(mjj) + 3I'{m H ) 2 ] 
A 2 

- — $ 2 [3K(m H , m H ,m H ) + (N - l)K(m H , m G , m G )) . (28) 

On dimensional reasons, and using the expected form of the divergencies we get for the relevant 
parts of the functions / and K (Don't forget that we deal with three-dimensional theory!): 

7(m) = 2J 1 Am 2 + 2J m 3 , 

A 2 2m + M 
K(m,m,M) = L log log — + L - 2Li og log . (29) 



Let us choose the counterterm 



5m\ = Ami ~ ~( N + 2) JiAA, (30) 
3 



what leads for the one loop part of (|28|) to the result: 



3/2 / \ \ 3 / 2 



vllL = \{™ 2 + Am 2 )$ 2 + A$4 + j Q p + + {N 1)Jq ^2 + j , , 

Using (pO) combined with the two loop counterterm 
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A 2 A 2 
dm 2 = Am 2 + — (N + 2)U og log — , 



(32) 



the two-loop effective potential (E8f) becomes: 
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2ren 
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A 2 



- ( m 2 + Am 2 + Am 2 - -(N + 2)L Q + -(N 2 + 8) J 2 ) $ 2 + 



A 
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3/2 



Jo m 2 + -$ 2 + (N - 1) J m 2 + 
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G 



3/2 



^Am 2 J ( m 2 + ^$ 2 



1/2 



A 



1/2 



A 



+-(iV - l)Am( Jo (^m 2 + -$ 2 j + _(JV - 1)AJ 2 ( m 2 + -$ 2 

^ 2 *? r To, 3m H 2to g + m/f 
+ ^ $ o^o S 3 log + (JV-l)log 

lo [I [I 



1/2 



m 2 + 



1/2 



(33) 



We can choose, in particular, the scheme 

Am 2 = 0, 



Am 2 = A 2 



iV + 2 
18 



Ln - 



N 2 + 8 
8 



Jn 



(34) 



which gives for the finite potential 



eff . 
2MS" 
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T ttt^o + Jo I ™ 



24 



3/2 



+ (N - 1) Jo ( m 2 + -$ 2 



3/2 



A 



+ _(7V_i) A J 2 (m 2 + -$ 2 
A 2 



1/2 



1/2 



3 log h (iV — l)log 



/./, 



/./, 



(35) 



Explicitly performing the 3D tadpole and setting sun calculations (see Appendix A), the values 
of the constants are: 



Ji 



4tT 2 ' 
1 



J °- 127T' 

Llog = 32^ 



6.70322- 10" 



(36) 



With these values 



eff _ ^2^2 , _^_<T> 4 L 
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m 2 + -$o + (Ar - 1) m 2 + -$ 2 
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3 log h (AT — 1) log 



(37) 



The divergence structure of the 3D theory can be read off from (30) and 
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(38) 



This coincides with the result of 0, if we change A — > A/6, which corresponds to the choice of the 
coefficient of the scalar selfinteraction to be A/4. 

The mismatch between the induced counterterms (E7p and the three-dimensional divergencies 



(38) is by now quite obvious. There are some divergences present in (27), which don't appear in 
(33) (eg. ~ ATlog(A/T)), but even the operators present in both cases have different coefficients 
(eg. ~ T 2 log(A/T)). This means, that in a regularisation scheme, where power-like divergences 
are simply subtracted (see ||), one still could recognize the need for supplementary terms in the 
action of the effective theory. The divergence structure of the reduced theory is richer anyhow 
than we could expect from the simplest (superrenormalisable) representation of the theory, eg. the 
scale dependence (beta- function) of the faithfully reduced theory is different. For the resolution 
of this puzzle, and for reconstructing the correct scale dependence we have to introduce some new 
operators into the action of the effective model. This will be proposed in the following section. 



IV. FILLING THE GAP: THE ONE LOOP FOUR POINT FUNCTION 



To cure the disease of the naively constructed 3D model detected in the previous section we have 
to add some new operators to the basic (||) action. The arguments in the Introduction suggest 
to investigate the four-point function (Q) and substitute it by a new vertex. Since, already all 
momentum-independent operators with couplings proportional to non-negative powers of T are 
taken into account, the new operators have to be momentum-dependent. From the point of view 
of the full theory, as we argued in the Introduction, the new term should represent the diagrams 
with mixed static and non-static internal lines. Such operators can be generated only radiatively, 
namely in the process of the integration over the nonstatic modes at 1-loop. In a subsequent 3D 

1- loop computation of its contribution to the effective potential, 2-loop divergencies yet missing in 
the comparison of the induced counterterms with the divergencies of the 3D local effective theory 
should emerge. 

We are going to evaluate the 1-loop correction to the 4-point function with static external lines 
of non-zero spatial momentum. When contracting two of its external lines as part of the solution 
of the effective model one recognizes the "mixed" static-nonstatic sunset contribution to the static 

2- point function. 
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The relevant integral (k = p\ + P2) is seen to be: 



(p 3 ,0)7f \ / f\(P4.,o) ,/A 1 1 

/A(fc) y p p 2 + (2irnT) 2 {p + k) 2 + (2irnT) 2 ' (39) 

where the primed integral is an abbreviation for the cut-off integration with respect to p 



(40) 



With Fcynman-parametrization one rewrites it as: 

1 

/A 



/A(fc) = ' dV I (p 2 + (2nnT) 2 + k 2 y(l-y)) 2 ^ (41) 



11 



9 



If the external lines are labelled by i, j, k ,1 respectively, the tensorial structure of the graph in 
Fig. 1. is proportional to 



(N + 4)%4; + 2(8^631 + S a S jk ). 



(42) 



Having this tensorial decomposition it is enough to work with homogeneous background as in the 
previous sections. 

The result of the integration is known for k 2 = 0: 



/a(* = 0) 



-Po-1 
8tt 2 



1 

8^2 



log — - log 2tt + 7_e - 1 



(43) 



In the other limit, k 2 3> T 2 , one can look for an expansion in positive powers of T (there is no 
other scale in the integral). Using the Euler-Maclaurin resummation formula: 



n r L 1 00 

£/(»)= dxf(x) + -[f(0) + f(n)]+J2 

i=0 n k=l 



E>2k 



[/ 



(2fc-l)( n ) _ j(2fc-l) 



'(0)] 



(44) 



(where i?2fc's are the Bernoulli-numbers) we get the following result: 



Ja(* » T) 
where 
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8tt 2 



F(x) 
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167T 



E 



(-ip(j-l)LB 2j /327r 2 T 



lllnll-- 



1. 



(45) 



(46) 



Terms proprotional 1/fc™ (n > 2) give constant contribution in three dimensions - these are the 
finite corrections to the four point function. The first three terms, however, might lead to diver- 
gencies when in a subsequent 1-loop calculation two of the external lines of Fig.l are contracted. 
If we want to separate only the divergent contribution to the effective potential upon contraction 
of two legs of Fig.l, we can interpolate between the relevant pieces of ( fl3| ) and p5| ) with help of 
the following expression: 



where C = 2ire 1 7E and 



fi(fc) = 



J A (fc) 



TT 2 Tk 

T 2 + k< 



1 

8^ 



5 log(l 



log^-n(fc) 



4C 2 T 2 



(47) 



(48) 



In section III, we have obtained for the one loop level momentum-independent radiative correc- 
tion (at k = 0) to the potential energy density from Fig.l: 



iV + 2 A 2 



315 



lo S ^ - log 2?r + 7 S - 1 



3 128tt 2 

This will now be modified for x-dependent fields $0 by adding the operators: 

N + 2 A 2 



AC 



3D 



1287T 



(49) 



(50) 



The Fourier-transform of ( p0[ ) promptly reproduces the k-dependent part of (f47j). For the modi- 
fication of the original theory (^) we have to project back ( ]50| ) the O(N) tensor-structure arising 
from (E||). The corresponding two possible O(N) symmetric operators 
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2 = E <Pitpjft(id)(pi<pj 



(51) 



have the relative weight N + 4 : 4, so for a general (/(-configuration we have to add the following 
operators to the basic action (||): 



A 2 f N + 4 4 
128^ — ° 1 + 9° 2 



(52) 



This is the central proposition of this paper. 

These operators are generated at one-loop level, so for two loop calculations of the effective 
potential they must be treated at one loop in the three dimensional theory. Whence the one loop 
contributions of different modes are independent we can calculate the extra divergences produced 
by these operators from the Lagrangian 



-3D 



N 1 



m%<$) + AC 



3D- 



(53) 



Shifting the fields tp i — > ipi(i ^ N) and (fx — » </?jv + 3>o respectively, and using the fact Q(k = 0) 
0, we get for the quadratic terms: 
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<%> = E \ (W 



i=l 
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2 



2 2 , 
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The effective potential using the above expression is: 



v eff = -mi$„ + (AT - 1)- / - — — In 
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2 J (2tt) 3 
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k 2 + m'i + 



(Af + 8)A 2 $ 2 
144vr 2 



n(k) 



After expanding it with respect to <& 2 the contribution to the quadratic term is: 

A 



,( 2 ) 
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"eff- 2 ' 



N + 2\ A 2 * 2 , 
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(54) 



(55) 
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The leading divergencies are provided by 



(2)dv 
"eff 



N + 2\ A 2 $ 



64tt 4 



dk 



, k\ ir 2 Tk (k 

log lr ) + t*T¥ + f a 



(57) 



what after performing the integrations simplifies to 
'N + 2\ A 2 * 2 , 



64tt 4 



A log A + 7T T log A - A In 47r + 1 - -y E 



, 2-i, 2 - a; 
aa; In — — 



(58) 
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The divergencies of the three dimensional model up to two loop level have to be cancelled against 
the induced counterterm Sm 2 nd (see (p7|)). There are divergencies coming from the local terms ( |38| ) 
and others coming from the nonlocal terms (|58|). Indead, summing these terms all divergencies, 
including the linear ones are cancelled exactly. This means, that for the proper treatment of the 
divergencies (or for the proper treatment of the scale dependence of the reduced theory) we must 
include the new operators (p2) into the effective description . These are nonlocal operators, the 
nonlocality is of the order OlT^ 1 ) in the coordinate space. On scales larger than C(T _1 ) these 
can be omitted. So a coarse grained local action with grain size C(T _1 ) may be a good choice to 
examine the whole theory R] - on this scale the theory forgets of its 4 dimensional origin. In the 
course of coming down from scale A to scale 0(T) in the momentum space the nonlocal operators 
modify the coupling constants of the local ones. The actual magnitude of this modification can 
be estimated only after doing this calculation, for instance, by matching some n-point functions. 
This is, however, the task of future calculations. 

V. CONCLUSIONS 

In this paper we advocated a non-local effective 3D theory, equivalent to the finite temperature 
O(N) symmetric scalar theory on the 2-loop level. The proposed theory has the form 

£eff = + m 2 {T)4> 2 a ] + + _^_[^±l^n(ia)^ + 1^0(0)^], (59) 

with m 2 (T) given by the coefficient of $o/2 in (p6f), and fl(id) defined in (ff8j). Now with a further 
step one can define a local approximation to this theory, when all static degrees of freedom are 
integrated out with k >> T. A local cut-off theory with A 3 = kT, k ~ 1 — 5 will be arrived at what 
is very close in spirit to the effective theory approach of [||. The simplest version of this theory 
has approximately the form 

£*fMf = \[Z{diKf + (M 2 (T) + E)#] + A3 ^ 4 A V , (60) 

where Z, E, A A3 should be found by matching some important quantities calculated from the 
respective theories ( ]59| ) and (|60|). We propose to calculate the couplings of ( |60[ ) by matching the 
low-fc behaviors of the 2-point functions and the self-coupling (4-point function). The matching is 
achieved the easiest way in the symmetric (high-T) regime, where m 2 {T) > 0. 

In the non-local theory our construction ensures the finiteness of the 2-point and 4-point func- 
tions. In the cut-off theory the dependence on the cut-off is not absorbed by anything, but is 
expected to be very weak if k is varied around unity. It will be very interesting to see the relation 
of the result of this approach to that of j(| . 

The practical interest of this second approximation is its easy discretisation into a lattice theory. 
However, since it is a cut-off theory, it should not be understood as a continuum theory, but it is 
better to study it with finite lattice constant of C(T _1 ) (coarse grained lattice). This means that 
the thermodynamics of this model should be studied with a dimensionless temperature = aT 
varying around unity. This circumstance seems to be very advantageous in view of the very weak 
phase transition signals, experienced in the small O region in some effective models ||. 

The analysis of the present paper should be carried over to gauged Higgs theories. The vector 
nature of the gauge potential does not present serious difficulty. More important is the task of 
analysing another class of diagrams contributing to the non-local behaviour of the effective theory, 
namely the four-leg box diagrams with non-vanishing external spatial momenta. 
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APPENDIX A: COMPUTING THREEFOLD CUT-OFF INTEGRALS IN 3 

DIMENSIONS 



We shortly discuss the method of performing 3D integrals relevant to the evaluation of the setting 
sun diagram contribution both in 3 and 4 dimensions. Here we use the regularisation convention 
that momenta on any internal line should not exceed the value of the cut-off, separately. (Another 
natural procedure would be to cut-off the loop momenta.) The integral to be performed is: 



d 3 p d 3 q d 3 k 
(2tt) 3 (2tt) 3 (2tt) 2 



6(A - |p|)0(A - M)6(A - \k\)(27r) 3 S^(p + q + k) f(p 2 , q 2 , k 2 ) (Al) 



Let us perform first the fc-integral with help of the delta-function - then instead of k 2 we can 
write (p + q) 2 . The remaining integral depends on p 2 , q 2 and x 2 — (pq) 2 /(p 2 q 2 ) (—1 < x < 1); 
k 2 = p 2 + q 2 + 2|p||g|a;. One can replace x by the variable k (its allowed range of variation is 
\p — q\ < k < p + q), then the integral will be symmetric in p, q, fc, with the domain of integration 



A:~{k<p + q p<q + k q<k+p}. 



(A2) 




Fig. 2 Integration range at fixed p 



Our expression after the angular integrations is of the form 



8^1 



dpdqdk 6(A - p)6(A - q)Q{A - k) pqk f(p A ,q A ,k 2 ). 



(A3) 



or writing out the limits of integration (cf. Fig. 2) 
1 



8tt 4 



\ A p+q A A p+q 

dp J dq J dk — j dp j dq J dk 

\ p - q \ A-p A 



pqkf(p 2 ,q 2 ,k 2 ) 



(A4) 



Often happens, that the integrand (/) is symmetric to the interchange of p and q. In this case we 
can write a simpler formula: 



/ = 



4tt 4 



A p p+q A p p+q 

dp J dq J dk — J dp J dq J dk 

p~q A/2 A-p A 



pqkf{p\q\k') 



(A5) 
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APPENDIX B: COMPUTATION OF THE SETTING SUN IN 4 DIMENSIONS 



The integral we have to evaluate is ( pi[ ) 

vt x f d 3 pd 3 qd 3 k 3 ^ ffi n+m+i (27r) 3 J(p + g + A:) mi x 

E(m p , m,q, m k ) = j^^T + ^ + ^ + E ^ (Bl) 

where E 2 = p 2 + m^. The summations can be performed with help of the " Saclay- method" Jl^] . 
We get the expression: 

K(m p ,m q ,m k ) = J ¥ * S(p + q + k)D{E p ,E q ,E k ), (B2) 



where 



2T 2 1 

D{E p , E q , E k ) = + AEpEqEk{Ep + Eq + Ek) 

T 



EpE q Ek 

1 



[Di (Ep, E q , Ek) + Di(E q , E p , E k ) + Di(E k , E q , E p )] 
D2[E p , E q , E k ) + D 2 (E q , E p , E k ) + D 2 (E k ,E q , E p )} 



and 



AE p E q E k 

- 7p I j-i [D 3 (E P , E q , E k ) + D 3 (E q ,E p , E k ) + D 3 (E k ,E q , E p )\ , 
/ i£jp£j q hj k 



D 1 (E p ,E q ,E k ) = ^( Tr ^—-—^—]+( q ^k) + 



(B3) 



2E p \E q + E k E q — E k J ' 2E p (E q + E k ) 

D 2 (E v ,E a ,E k ) = n E | - - ), 

ZV p, K) Fj v I 771 I 771 I 77? 771 77 771 / 1 

\£jp + Ej q + tj k £jp — £j q — £j k J 



D 3 (E p , E q ,E k ) = n Eq n El 



p 1 ^q \ J-Jfc *-Jp ^q 

1 1 1 



E p + Eg + E k Ep — E q + E k E p + E q — E k 
1 



E p — E q — E k 



(B4) 



and tie is the Bose-Einsein function (tie = (e@ E — l) -1 )- The parts depending on T only in the 
form n e agree with the expression of Parwani Jl5| . 

In the course of the integration of each single term in (B2) UV and various IR divergence 
problems arise. In order to control them the best is to work from the start with a regularized 
propagator in the mixed (r, p) representation fl4|| : 

A reg (r, E p ) = A(r, E p + e)0(A - \p\), (B5) 

where 

A(r,X) = ||(e^ + e ^). (B6) 

With this propagator, using the general scheme to perform threefold integrals in 3 dimensions 
(Appendix A) we get the following expression for (|B2|): 

K reg = Jdpdqdk pqk 9(A - p)Q(A - q)Q(A - k)D(E p + e, E q + e, E k + e). (B7) 
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The full expression is IR safe, since we omitted all static parts. Therefore at the end of the 
calculations we may perform the e — > limit. 

If we rescale all the variables by T, the full integral K ~ T 2 . At high temperatures it is possible 
and convenient to expand the integral with respect to the masses: 



d 

K reg (m p , m q ,m k ) = T 2 K reg (0, 0, 0) + m 2 p -^—^K reg {m p , m qi m k ) 



+ (p «-> q) + (p <->■ k) 

m ; =0 



(B8) 



The omitted parts are proportional to m 4 /T 2 . 

The regularization of the propagator regularizes all the arising divergencies, which makes pos- 
sible the correct derivation of the finite parts of the integrals. This is important in the case of 
K reg (Q, 0, 0), where the finite part is multiplied by T 2 , which means, that it cannot be absorbed 
by counterterms in any T-independent renormalization scheme. 

For the calculation of K reg (0, 0, 0) I followed the following general ideas. First one has to 
separate the (UV/IR) divergent parts. They should have the simplest possible form allowing 
anlytic evaluation of the coefficients of the divergencies. Substracting the divergent pieces we can 
perform the A — > oo and the e — > limits respectively. If there is some Bose-Einsein functions 
(tie) involved in the integral containing A or e, we can use the relations: 



lim 

£^0 



lim 



dxen(ex) f{x) = 
dxA 2 n Ax f{x) 



dx 



fix) 



6 



f(0) 



if it exists, 
if/(0) = 0. 



(B9) 



After this step exclusively integrals over rational functions remain to be performed. The k- 
integration can be done analytically. The p,q integrations have been transformed to a form (in 
general by introducing new integrational variables), where at most 1 numerical integration was 
left. 

Armed with these ideas all the integrals in K reg (0, 0, 0) can be performed. Let us denote the 
integrational measure in (|B8| ) after shifting all the variables by e by 



p+q-e 



J dM c = Jdpip ~ £ ) fdqiq - e) J < 



p+q 



dp{p — e) J dq(q — e) J dk(k — e) — 

|p— flfl+e 



dpp J dqq 

\-p 



dkk 



(B10) 



In the second part of the measure no IR divergencies occure, the e — > limit can be taken directly. 
The results of the relevant integrals are (ignoring terms, which vanish as A — > oo or e — » 0): 



I A 

Ib 
Ic 
Id 
Ie 
If 
Ig 



dM e 
dM £ 
dM E 
dM e 
dM £ 
dM e 
dM F 



1 



1 



p 2 q 2 k 2 
1 



16tt 2 



In - - 0.0190158229, 



2pqk 2 (p + q) 



= 0.0009950146 A, 



2pqk 2 



p 



p + q 



p-q 



1 



16-7T 



In- -0.0115413434, 



1 



4pqk(p + q 

Up 

Apqk{p + q - 

n p 

4pqk(q + k — p) 
1 



k) 



k) 



0.0001041333 A 2 , 
1 



1927T 2 

1 

1927T 2 

1 



In A- 



0.0000776673, 

- 0.0003657838, 
1 



Apqk \k + p + q k + p — q k — p + q 



p-q 
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64tt 2 



■ mi 



0.0044038354. 



(Bll) 



With these integrals we can express we can express the first term in the expansion of (B8): 

K reg (0, 0, 0) = 21 A - 3I B - 3I C +Id + 3J b + 3I F + 37 G (B12) 

The IR divergencies cancel each other as we expected, the rest gives the result appearing in the 
expression (|25|): 

A 



^reg (0,0,0) = 0.0001041333 A - 0.0029850437 AT + — — log - - 0.0152887686 T 



32tt 2 



T 



(B13) 



Now, it is worthwhile to give a few commnets on the second term of the mass-expansion. Its finite 
part can be absorbed into the finite renormalization, so only the divergent pieces are important. 
Therefore less terms are to be preserved. The differentiation with respect to the mass squared can 
be changed to differentiation with respect to the energy or (since we put afterwards all masses = 0) 
to the momentum under the integral: 
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D 



A,e 



1 3D 



A,e 



2Ek dEk 



(B14) 



When the derivation acts on n(E), we can integrate partially, so at the end only rational functions 
and eventually n(E) factors remain, and with the previous ideas the integrations can be reduced 
to one-dimensional numerically computable integrals. 
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